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Abstract 

We obtain a general rule that the 0(1/ log m ./J term due to the current hy^y^h 
of a massquark h in /-flavour theory is —3g 1 2 3 4 (m/i)/{27r 2 (33 — 2/)} times 
the flavour singlet current of the residual (/ —l)-flavour theory, where gj is the 
/-flavour running coupling constant in a mass-independent renormalization 
scheme. The rule is applied to the Ellis-Jaffe moment below and well above 
charm threshold, and to low-energy Z°-exchange amplitudes. The singlet 
axial charge of the proton common to these experiments is both scale and 
gauge invariant, but is related to the axial anomaly and the “gluon spin” by 
a non-perturbative renormalization factor. 
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With major experiments underway to determine the strange-quark contributions 
to polarised deep inelastic scattering and low-energy parity-violating elastic ep and 
up scattering, it is vital to understand precisely what will be learnt. At first sight, 
the answer seems clear: throw away contributions from heavy quarks, including 
charm, and identify the remaining flavour singlet and non-singlet contributions to 
axial charges of the proton in the various reactions. 

However the situation in QCD is far from obvious because of renormalization 
effects in the singlet channel [|T], |2|. 3 associated with the axial anomaly 0, ||. In 
separate studies of the deep inelastic 0 and elastic neutral current 0 processes 
mentioned above, it has been noted that these effects give rise to a characteristic 
factor[] 

E(g) — exp / dg'p(g’)/(3(g') (1) 

Jo 

multiplying the gauge-invariantly renormalized singlet axial-vector operator 

J S = («7m75 u + dy M y 5 d + sy^s) (2) 

in the residual theory (only three flavours u, d, s if c, b , t are considered heavy). Here 
/3(g) and y (g) are the Callan-Symanzik functions of the residual theory associated 
with the gluon coupling constant g and the composite operator J^if. We note that 
the factor © is non-perturbative: it involves integration over a finite range from 0 
to a value g defined at a renormalization point /j well below the charm threshold. 

Our aim is to pin down the relation between the axial proton charges which arise 
in these inelastic and elastic processes, while taking full account of complications 
such as the factor (!l|). This turns out to be possible for leading terms in the heavy- 
quark expansion 

rrit, m b , m c —> oo , > m/, > m c (3) 

and includes a simple rule for the 0(1/ log rrih) correction as each quark h = t, b, or 
c is made infinitely heavy. For the inelastic process, two limits are considered: 

(a) Asymptotia below charm threshold, defined by 

rn t m b » m c Q 3> p (4) 

where /j is the renormalization subtraction point, i.e. all heavy quark limits, 
including m c —*■ oo, are taken before the limit Q = \J—q 2 —> oo of large 
space-like momentum transfer cp,. 

(b) At the end of the paper, asymptotia well above charm threshold, but below 
bottom threshold, viz. 


m t » m b » Q » m c » p (5) 

5 Many earlier papers, from jl|] or || onwards, make implicit reference to this factor but then 
either omit it or keep only the leading order in the exponent. The lower limit of integration in 
equation (9) of § is missing; presumably 0 was intended. We follow the convention of § for 7 (g), 
i.e. the two-loop term is positive. 
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We begin with inelastic polarized ep scattering in the limit (J|), the residual 
theory being QCD with three flavours u,d,s. According to (£|), asymptotia for a 
structure function W is the deep inelastic limit Q —> oo taken after the limit (g), 


W (. x , oo) 


lim lim W(x, Q 2 ) 

Q —>ooeq (3) 


6 flavours 


( 6 ) 


where x = Q 2 /2p.q is Bjorken’s variable. This means that the relevant deep inelastic 
data from which a three-flavour singlet charge is to be extracted should be restricted 
to kinematic regions below the charm threshold (A C D production): 

Q 2 (l-x)/x < {m D + m Ac ) 2 ~ m 2 p (7) 

Small violations of this condition may be acceptable, provided that all charm thresh¬ 
old effects are isolated [jSJ and removed. 

Let g = g(Q) be the effective coupling constant of the three-flavour theory, 


l°g (Q/v) = I dg'/P(g') 

J 9 

and let g\ | inv be the axial charge of the operator [6| : |;j 

(x) = E(g)J^/(x) 

for proton states | p, s) with momentum p p and spin s 

(p, s|S„(0)|p, s) = 2m p s il g° A . 


( 8 ) 

(9) 


( 10 ) 


Then, as Q tends to oo, the SU(3) flavour singlet part of the Ellis-Jaffe moment 
ig is given by g |. g 0, 0 


dx gi(x, Q 2 ) 


singlet 


1 .0 
9 


~ -,9a 


{i + Ec,9 2 '(Q)} 


( 11 ) 


i>\ 


where the coefficient q is calculable in Aloop perturbation theory g . 

The factor E(g) in (fl) arises naturally from the coefficient function of Jjff in the 
product of electromagnetic currents J a (x)Jp( 0) at short distances ~ 0. It com¬ 
pensates for the scale dependence of .J p f caused by the anomaly in its divergence^, 

f 

= 2 fd^Kp + E 2^*75® (12) 

i= 1 

where Kf is a renormalized version of the gluonic Chern-Simons current, and the 
number of flavours / is 3. 

6 As Adler noted in his original paper fl|, the anomaly gives rise to a multiplicative factor 
Zrj which makes jfff finite at two-loop level and beyond. The corresponding 7 (g) function was 
introduced by Koberle and Nielsen |l4j in quantum electrodynamics (before asymptotic freedom). 
The extension to QCD was discussed in the review |0 and independently in g. There followed 
Kodaira’s well-known application || to polarised deep inelastic scattering. 
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A consequence of the renormalization-scale invariance of is that its spatial 
components have operator charges 


Si(t) = Jd 3 xS l (t,x) 


which satisfy an equal-time algebra 


[Si(t),Sj(t)] = i£i jk S k (t ) 


(13) 


(14) 


characteristic of spin operators. Unfortunately, the time dependence of the “spin” 
Si implies that energy eigenstates such as | p, s) are complicated mixtures of its 
representations. Note that a proposal [|16|] to obtain the U( 6 ) x U( 6 ) algebra |L7| in 
QCD works only if the singlet axial-vector operator is identified as <S /( . 

As noted above, the factor ([I]) is not perturbative, and so should not be ap¬ 
proximated by retaining just the lowest-order terms in /3 and 7 . We also note that, 
although the divergence of is 

(15) 


= 6 E(g)d ll K^ + mass terms 


because K fl is not 


the quantity E(g)d ti K fi is not scale invariant (contrary to IB 
multiplicatively renormalized^: its counterterms are all proportional to ,/W. As a 
result, in the relation between the observed singlet charge g A | inv , the “gluon spin” 
Ag |L 8 ], |l]| as interpreted in the light-front gauge [|2(J, and the scale invariant 


contribution AE cons from the gauge dependent, partially conserved chiral current, 
E{g) appears as an uncontrolled theoretical factor:[] 


9a 


E(g){ AS cons - Ag) 


8n 2 


(16) 


The scale dependence of Ag shows that the canonical basis for supposing it to be 
gluonic spin is destroyed by renormalization. We emphasize that the results of this 
paper do not depend on assumptions about E(g) or Ag. 

So, returning to ([TT]) , we arrive at a conclusion which is now essentially con¬ 
ventional: the three-flavour singlet Ellis-Jaffe moment (|TT|) at Q ~ 00 tends to a 
singlet axial charge 7 ^ | inv which is both gauge and scale invariant, as any measurable 


quantity should be. Its current experimental value is [21, 22 


9a 


= 0.28 ±0.07 


(17) 


Other notations for this quantity are E inv || and AS( 00 ) 

Now let us consider whether 0 this deep-inelastic version of the singlet axial 
charge is relevant for elastic Z°-exchange reactions. Once again, we will consider 
the heavy-quark limit ( 0 ), in this case, applied to the axial-vector neutral current 


T z - 

J fi5 — 


E - E 

q=u,c,t q=d,s,b 




(18) 


'This point has been rediscovered several times. Again, the essence of it is in Adler’s original 
paper [Q, being connected with non-renormalization of the anomalous constant. The minor changes 
needed for QCD were discussed in the review & See the reviews jTT[ [Ui| for recent discussions of 
this point in the context of proton “spin”. 

8 This reflects the well-known fact that the ratio of the moment (0) at Q = 00 to that at low 
Q is non-perturbative. Our point here is that we will be able to relate deep inelastic scattering to 
elastic Z° exchange despite the presence of such a factor. 
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coupled to the proton. 

Originally, Collins et al. [|lj estimated residual effects from anomaly-cancelling 
pairs (t,b), (c, s), ...in ([18|) . Later Kaplan and Manohar [^] considered how to 
extend the method to the more realistic case where c is heavy and s is light. The 
procedure adopted in these papers was to set a renormalization scale /q n t between 
each quark made heavy and the next in line. This approach belongs to the class 
of mass-dependent renormalization schemes in which decoupling theorems |25| are 
most easily proven. Comparisons with mass-independent schemes, where decoupling 
is not manifest, were the subject of much discussion many years ago p 6 |, pT| . 

However, mass-independent schemes are preferable when one wants to develop 
explicit asymptotic expansions in 1/m/,.. For example, Chetyrkin and Kuhn [[7], 
commenting on the result of JIJ, observed that the renormalization prescription 
of amplitudes in the residual theory produced by a heavy-quark limit is hard to 
determine in a mass-dependent renormalization scheme. It is much simpler to have 
asymptotic results depending on just one scale parameter /i and not on the residual 
quark masses. 

So we choose to work in a mass-independent scheme, choosing /r such that it 
is suitable for the final residual theory with just three flavours, and relying on the 
renormalization group to keep track of all factors. This means that we consider /i 
to be well below m c , as in the case of deep inelastic scattering. 

Let Aq be an amplitude in the full six-flavour theory. Its expansion in the inverse 
top mass rrii * looks like an operator-product expansion in momentum space, but 
with Q replaced by m t . The result involves six-flavour coefficient functions C n e(m t ) 
which multiply five-flavour amplitudes A n 5 : 


A 6 ~ ^ 2 C n 6 (m t )A n5 , mt —► 00 

n 


(19) 


Since all C n d(m t ) are obtained at rn t ~ 00 , they correspond to sub-amplitudes in 
which all light-quark (q 7 ^ t) and gluon propagators carry loop momenta 0{m t ). 
This means that, provided the renormalization scheme is mass-independent, these 
coefficient functions do not depend on mb or lighter quark masses. 

Consequently, the heavy -6 expansion can be obtained by expanding each five- 
flavour amplitude A n 5 in ([L9D in terms of five-flavour coefficient functions C nn /5 and 
four-flavour amplitudes A nn / 4 , 

And {jTlb)Ann'A t ( 20 ) 

n' 


and substituting back into (jT9[) . The limit sequence can be continued by then 


expanding A nn 'i in m' 1 , but no further, since c is the lightest quark which can 
conceivably be treated as heavy. 

The full expansion of the original amplitude (pT|) is thus 


Ad 


) ) Cnd (’^t)C nn '5 (nib)Cnn l n"4 ij^c)A n n'n"‘i i Wit Wlb /i (21) 


It involves non-perturbative three-flavour amplitudes A nn ' n "3 which we would like 
to relate to quantities such as g^| inv . 


4 




Let the six-flavour amplitude (® be the proton matrix element of the neutral 


current ([T8|) : 

A = ip,s\J^\p,s) ( 22 ) 

It is well known that J^ 5 has a soft divergence and so is renormalization-scale invari¬ 
ant in the full six-flavour theory. Consequently, the amplitude (|22|) is scale invariant, 
as is required for a measurable quantity. For any heavy-quark limit, it must remain 
measurable and scale invariant in the residual theory with fewer flavours. 

First we consider terms in (j2l|) which survive the limiting procedure; logarith¬ 
mic corrections will be considered later. The residual theory has three flavours, so 
flavour SU(3) should be used to classify operators in the expansion of J^ 5 . The only 
possibilities are the octet operators 


= 9^7sA a g 


(23) 


and the singlet operator (| 2 |), or better, the scale invariant operator S defined in ([]), 
since we know that the answer must be scale invariant. Clearly, the terms q = u,d in 
(|TB|) give J^ 5 , the term q = s can be split, into — J® 5 /a/ 3 and something proportional 
to Sfj , while the terms q — t,b,c contribute to only: 

J f f 5 —> octet + {constant}^ (24) 

By construction, the constant of proportionality in ( p4j ) is scale invariant. It is also 
quark-mass independent, since we use a mass-independent renormalization scheme, 
so it must be independent of the gluon coupling constant g. Therefore we can 
evaluate it by considering the limit g —> 0 for J- f 5 coupled to off-shell u, d, s quarks 
and noting that F'(O) = 1: 

constant = — - (25) 

6 

We conclude that, in the strict heavy-quark limit (|3]), the neutral current is given 
by 

4 —> - \S r (26) 

Consequently, there is a precise relation between the three-flavour Ellis-Jaffe 
moment ([II]) and the neutral current amplitude (^) in the limit (Jl), 


{P,s\J%i\p, s ) —* 2 rripSf, (gf + -^gf - ^g° A 

where g^ are the usual octet axial charges 

(p, sl-'JIP. s ) = 2m, a = 3,8 


(27) 


(28) 


(so gA/gv for neutron decay equals 2 g\). Note that (|27| ) holds despite the presence 
of E(g). The matching condition (^5|) at g = 0 is crucial in arriving at this otherwise 
conventional conclusion. It can be made to look even more conventional, 


{p,s\J. 


\Pi S ) 


(2m 


pS/i, 


2 (^A^inv d[ nv As, 


(29) 
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if we define scale invariant axial charges Ag inv for q = u,d, s via the prescription 


Alt i nv 

1 

CO 

0 . 

CO 

( & ) 


Ad; nv 

— 3 

-3 73 1 

9a 

(30) 

/ A Sinv j 


^ 0 -2^3 1 ) 

V sSU 1 



Apart from its strong dependence on the heavy-quark limit, the relation between 
(|TI|) and (p7| ) is an isoscalar version of Bjorken’s sum rule f24jl , 

The next step is to determine the leading 0(1/ log rrih) corrections to (p7|). 

Let us retain only the leading power in the general expansion ([HI) . If we arrange 
that each independent five-flavour amplitude A$ remaining on the right-hand side 
is scale invariant, the result necessarily consists of one or more expressions of the 
form 

A.6 ~ ^2 Cm + • • • (31) 

£>0 

where ge is the effective gluon coupling in the six-flavour theory and the coefficients 
Cm can be calculated perturbatively. Proceeding as before, we find that the leading 
power of (pID consists of one or more terms 

Aa ~ A 3 J2 CmC Vb Cmigl\m t )gf'(m b )ff(m c ) + ... (32) 

e,e,e">o 


involving the effective gluon couplings gy for / flavours as functions of rn t ,m b ,m c 
for / = 6,5,4 respectively. (In this notation, the effective coupling g in ([□]) would 
be written g 3 (Q).) 

Suppose that, at a given stage of the calculation, / flavours remain, and h is the 
next quark to be made heavy. As noted by Collins et al. [| 1 ] , the 0(1/ logm/,) correc¬ 
tion arises because the current h'y^^h in ,/// can be part of a triangle sub-diagram 
formed from h-quark propagators. For rrih large, hj^'y^h mixes with the operator 
Hqjkh 177 / 05 Q- We need the scale-invariant versions of these currents contained in the 
column vector 


V = ( HOs h \ 

M V 97/059 ) ■ 

where the equations 


r 1 -1 
1-r 1 1 


{S»}f \ 
{^ 5 }/) 


{S,}f = {E(g)}f{i: GJ 

2—1 

/ 

i Y ^}f = l" 1 XoO/05^ - HOs^ 

2—1 


(33) 


(34) 

(35) 


dehne /-flavour analogues of ([|) and the axial-vector hypercharge current. Note 
that the currents in V M are normalized such that at g = 0 , they reduce to the usual 
free-held expressions. 

Let 7 (g) be the matrix Callan-Symanzik function controlling the mixing of 
components as rrih becomes large. At rrih = 00 , the lower component of tends 


6 








to which is the analogue of (34) in (/ — l)-flavour theory. Therefore the 

entire leading powei^j is given by 





{hiffoh} inv 


+ 0(m 


-i\ 


m,h =oo 


(36) 


where V denotes ordering of the exponential. 

Our interest is in the upper off-diagonal term g 4 / 167 r 4 + 0(g 6 ) in 7 combined 
with —g 3 (33 — 2 /)/ 487 r 2 + 0(g 5 ) for the /-flavour (3 function (for three colours). 
Substituting these expressions in (|3b|), we deduce our most interesting result, viz. 
the rule 

{^ 75 ^}inv = ~ 2^33^/) + ( 37 ) 

for the leading 0(1/ log m/J correction. This analytic result is the hrst of a series of 
perturbative calculations which can now be done to determine various coefficients 
in (|32|). The techniques resemble those used to obtain 0(1/log Q) expansions such 
as flU]). 

We should note that Chetyrkin and Kuhn [^] come very close to stating this 
rule, but they differ on a detail in the flavour dependence. This becomes clear 
when successive heavy-quark limits are considered. Following [j|l|, they considered 
m t —> 00 and then mj —> 00 for the difference i'y fl 'y 5 t—b'y fl 'y 5 b and found a /we-flavour 
expression 

|(j!W 

for the coefficient of the four-flavour singlet current. Our rule requires a six-flavour 
dependence for the 0(1/ log'/rq) term, / = 5 being relevant only for the 0(1/ log up,) 
term: 

7 7m75 t - bvlsb ~ (Jj^f (™ft) - \gl(m t )) ^2 {^} 4 (38) 

ft is not clear whether the result of Collins et al. [f] is consistent with ( PB| ) or 
not. Their formula contains only one running coupling g evaluated at both 777 and 
mb, so the flavour dependence does not appear to be the same. However, their 
renormalization scheme is mass dependent (unlike ours), producing amplitudes in 
the residual theory with renormalization prescriptions which are hard to specify [0. 
Comparisons with the work of Kaplan and Manohar [;2|. |28|, in which the mass- 
dependent approach of Q is taken much further, are correspondingly more difficult 
to make. 

In the limit @, the rule ([?7]) applied to the neutral current implies 

~ + 0(l/m„) (39) 


where the coefficient C (including the factor 


in 


(HD) 


IS 


C = \ + ^2 + ^gl(m c )) + 0({log(m fc )} 2 ) (40) 

9 Strictly speaking, we should replace the renormalized parameter mh by the scale invariant 
quantity Km^fexp — // ^mig')/ Pig 1 )} f in all asymptotic equations, where 7 m {g) is the Callan- 
Symanzik function for mass renormalization in the /-flavour theory, and k is a number to be 
determined from the finite part of the two-loop calculation at mh = 00. This does not affect 
0(1/ log mh) corrections. 
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This coefficient fixes the first logarithmic correction for the amplitude (P7|): 


(p,s \ J Zs\p,s) = 2 m p s^ 


gf + -^9 { a ~ c ( m h) g°A 


+ 0(l/m h ) (41) 


As in (g5j), normalizations can be checked via matching conditions at g ~ 0 at any 
point in the derivation. 

A numerical estimate of the corrections in C can be obtained by substituting 




247r 2 


33-2/ 


{log(m h /n)} 


-i 


(42) 


with the scale g chosen well below m c but in the range of a few to several hundred 
MeV characteristic of low-energy strong interactions. Let us take g = 0.5 GeV. 
Then, listing contributions in the order h = t, b, c, we find 


C « i(l+ 0.04 -0.09 + 0.15) = 1.11) 


(43) 


The correction is almost three times that of @ and opposite in sign. It corrects As inv 
by an amount 0.01 compared with the experimental value |U], —0.11 ± 0.03 for 

inelastic processes. 

The rule ([U]) can also be used to compare the Ellis-Jaffe moment (|TT]) below 
charm threshold, as defined by the limiting procedure (|4[), and the corresponding 
quantity Jq dxg^(x, Q 2 ) well above threshold, where asymptotia is defined by the 
condition (|^). Note that c is still considered heavy, not light (since flavour SU(4:) 
symmetry is a very bad approximation compared with m c g). Going from ([|) 
to ([5]), we have to replace the three-flavour singlet operator S by the four-flavour 
expression 

{5Z gi^g + 2 c7 At 75c}. nv (44) 


q^c 


Consequently, as m c becomes large, the rule (0) predicts 


J dx (gi( x , °°) - gi( x , oo)) — T^glimc) 


la a 


h a 


(-0.10) (45) 


where S'Aliny is the three-flavour singlet axial charge ([[(]) appearing in ([□]) and (|ll|). 
The minus sign in ([45| ) indicates that charm production is slightly enhanced when 
the beam and target spins are aligned. 

A previous estimate [^] of (4S) was obtained by evolving the moment through 
the charm threshold. This involves a specific model for threshold behaviour for 
scale-dependent amplitudes in a mass-dependent renormalization prescription, so a 
theoretical or numerical comparison with our result is not obvious. Throughout, 
we avoid thresholds and use a mass-independent prescription; that leads to explicit 
results for 0(1/logm/J corrections. 
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